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Abstract: 

We present a fit to spin asymmetries which gives polarized quark distribu¬ 
tions. These functions are closely related to the ones given by the newest 
Martin, Roberts and Stirling fit for unpolarized structure functions. The 
integrals of polarized distributions are discussed and compared with the cor¬ 
responding quantities obtained from neutron and hyperon /3-decay data. We 
use the combination of proton and neutron spin asymmetries in order to de¬ 
termine the coefficients of our polarized quark distributions. Our fit shows 
that phenomenologically there is no need for taking gluonic degrees of free¬ 
dom into account. 
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The new data taken in experiments at CERN |I. 2j] and SLAC |3j] renewed 
the interest in the spin structure of the nucleon. We have now at our dis¬ 
posal relatively precise CERN data for proton spin asymmetries from Spin 
Muon Collaboration (SMC) || and very recent E143 experiment results from 
SLAC [ff|. We have also the asymmetries measured from 3 He in E142 exper¬ 
iment at SLAC [Q whereas the data from deuterium target come from SMC 
0 ], Together with an old SLAC || and EMC |§ data for proton one has 
considerable amount of information which can be used to study nucleon spin 
structure and in particular to determine the polarized quark distributions. 

The unpolarized quark distributions in the nucleon are known due to 
several fits [7-10]. Martin, Roberts and Stirling (MRS) @ gave a complete 
fit using an existing experimental data in order to determine parton (i.e. 
quark and gluon) distributions. Quite recently an improved version of such 
densities has been presented by them ||. 

Our preliminary discussion how to determine the polarized quark distri¬ 
butions using the older version D_ of the MRS fit was given in ||. 

In this paper we would like to get polarized quark parton distributions 
starting from the unpolarized ones and using existing data for proton, neu¬ 
tron (made on 3 He) and deuteron spin asymmetries. We will consider two 
different models of possible x behaviour (more and less singular at x —> 0) of 
the sea contribution. The calculated values of octet axial-vector couplings a 3 
and a§ obtained from the fit are compared with the experimental values got¬ 
ten from nucleon and hyperon /3-decays modified for QCD corrections. These 
quantities are used to differentiate between various fits (together with their 
X 2 values) which were made for different combinations of the spin asymmetry 
data (e.g. proton-)-neutron versus proton+deuteron data). Putting the addi¬ 
tional restriction for as (in order to stabilize the fits) we determine polarized 
quark distributions from the combination of proton and neutron data. We 
have also tried to include polarized gluons contributing in the way proposed 
in Ref. [IT]. It does not lead to any substantial improvement (y 2 per degree 
of freedom is worse) in the fit. In addition the sign of the gluon contribu¬ 
tion is opposite to the one which is expected theoretically. The inclusion of 
newest SLAC data for proton asymmetry from the E143 experiment is also 
discussed. 

Let us start with the formulas for unpolarized quark parton distributions 
(at Q 2 = 4 GeV 2 ) given by Martin, Roberts and Stirling ||. We shall consider 
their fit called MRS (A) with rather singular behavior of sea distribution at 
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small x values (which, however agrees with the results from HERA 
have for the valence quarks distributions: 

u v (x) = 1.996x”°- 462 (1 -x) 3 ' 96 (1 - 0.39^ + 5.13x), 
d v (x) = 0.296x- a670 (l »x) 4 - 71 (1 + 5.03Vx + 5.56x), 

and for the sea ones: 

2 u(x) = 0.392 M(x)-S(x), 

2 d(x) = 0.392 M(x) + S(x), 

2 s(x) = 0.196 M(x), 

2 c(x) = 0.020M(x). 

In eq.(2) one has for singlet: 

M(x) = 0.411oT L3 (l - x) 9 ' 27 (l - 1.15y/x + 15.6a:), 

and isovector part: 

S(x) = 0.099 x~ OA 62 (1 - x) 9 ' 27 (l + 25.Ox). 

The unpolarized gluon distribution is given by: 

G(x) = 0.775x _ 1 ' 3 (1 - x) 5 - 3 (l + 5.2x). 


We 
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In analogy to the unpolarized case we assume that the polarized quark 
distributions are of the form: x"(l — x) /3 P 2 (y/x), where P 2 (\/x) is a second 
order polynomial in yfx and the asymptotic behavior for x—>0 and x^l (i.e. 
the values of a and (5) are the same as in unpolarized case. The unpolarized 
parton distribution is the sum of helicity distributions along and opposite 
to parent nucleon helicity whereas the polarized one is the difference of such 
functions. Hence, our idea is just to split the numerical constants (coefficients 
of P 2 polynomial) in eqs. (1,3,4) in two parts in such a manner that the 
distributions are positive defined. Our expressions for A q(x) = q + (x)—q~(x) 
(q(x) = q + (x) + q~(x)) are: 

A u v (x) = a;^ a462 (i _ x) 3 - 96 (ai + a 2 \fx T a 3 x), 

A d v {x) = a;- 0 - 670 (i - x ) 4 ' n (bi + b 2 \fx + b 3 x), 

A M(x) = x~ 0800 (l — x) 9 ' 27 (ci + c 2 y/x), (6) 

Ah(x) = x^°' 462 (l — x) 9 ' 2 'd(l + 25.0x). 
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At this stage we will not take into account gluonic degrees of freedom. In 
the case of AM, i.e. total sea polarization, we assume that there is no term 
behaving like £~ L3 at small x (we assume that AM and hence all distributions 
are integrable), which means that coefficient in this case have to be splitted 
into equal parts in M + and M~. The next term (;r~ 0 ' 8 ) is relatively singular 
in comparison to valence quark distributions. That means that for x —> 0 the 
sea contribution dominates and hence, proton and neutron spin asymmetries 
will behave in the same way in this regime. That is not the trend that is 
observed in the existing experimental data. Looking at the data points we 
see that proton asymmetry is positive while neutron negative for small x 
values. Nevertheless, we shall consider two kinds of models: the first (later 
abbreviated by I) with more singular (AM ~ x~ 0 ' 8 ) behaviour for x — > 0 
and the second (II) with less singular (AM ~ ;r~ 0 ' 3 ). In the second case the 
coefficient in front of a; -0 ’ 8 is also equally divided between M + (x) and M~(x). 
We would like to stress that such behaviour is not a Regge type extrapolation 
in this region. There are some problems in MRS (A) fit with positivity of 
quark sea distributions and the condition that q + and q~ distributions have 
the same minimal non positive term determines d in eq.(6) and reduces the 
number of parameters from 8 to 7 (or from 9 to 8 when c± % 0). 

To find the unknown parameters in the expressions for polarized quark 
distributions we fit our formulas for Ai, A™ and Af (with eight or seven 
parameters) to the experimental data on spin asymmetries, which are given 
by: 


A\(x) 

A^x) 

Af(x) 


4A u v {x) + A d v (x) + 2.236A M(x) — 3A5(a;) 

4 u v (x) + d v {x) + 2.236 M(x) — 35(x) 

A u v (x) + 4A d v (x) + 2.236A M(x) + 3A 5(x) 
u v (x) + 4 d v (x) + 2.236M (x) + 35 (x) 
5Au v (x)+4.472AM (x) 3 

5n v (x) + 4.472M(x) (1 “ 2 Pd){1 + R) ' 


(1 + R), 
(1 + R), 


(7) 


The ratio R = ctl/ct, which vanishes in the Bjorken limit, is taken from 
Rj, whereas po is a probability of D-state in deuteron wave function (equal 
to 5.8%). Spin structure function g\ is given by: 


gi(x,Q 2 ) = (4A-u„(o;) + A d v (x) + 2.236AM(x) — 3A<5(a:))/18. (8) 


In this paper we assume that the spin asymmetries do not depend on 
Q 2 what is suggested by the experimental data IB and phenomenological 
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analysis |L4|]. We have made fits using different assumptions about the sea 
contributions (I and II) and also for different sets of experimental data (i.e.: 
proton+neutron (pn), proton+deuteron (pd) and proton+neutron+deuteron 
(pnd)). We do this because it is known that neutron and deuteron spin 
asymmetries lead to different values of an integral: T”(Q 2 ) = Jq g™ (x, Q 2 ) dx. 
We have also tried to include gluons along the line of ref.[|ll|] by considering 
effective polarized quark distributions A q e ^ = A q — |^A G and fitting the 
constant in front of gluonic distribution. 

The obtained polarized quark distributions A u(x), A d{x), A M(x) and 
A<5(x) can be used to calculate hrst moments. For a given Q 2 we can write 
the relations: 


r? = — Au + — Ad+ — As, 

1 18 18 18 

14 1 

r? = —A u H- Ad H-As, (9) 

1 18 18 18 ’ y J 

where A q = /q A q(x) dx. 

Other combinations (octet and singlet ones) of hrst moments of quark 
polarizations are: 


03 = A u — Ad, 

a 8 = A u + Ad — 2As, (10) 

AS = A u + Ad + As, 


The results for the integrated quantities (calculated at 4 GeV 2 ) after tak¬ 
ing into account known QCD corrections (see e.g. Ref. |[HJ ) could be com¬ 
pared with axial-vector coupling constant Qa and g$ known from neutron 
/3 -decay and hyper on /1-decays (the last one with the help of SU(3) symme¬ 
try). We can express the combination of T^Q 2 ) and T"(Q 2 ): 

r?(Q 2 ) - r?(Q 2 ) = c NS (Q 2 )g A /6, (11) 


where cns{Q 2 ) describes QCD corrections for non-singlet quantities [T5j and 
g A = 1.2573 ± 0.0028 (see Ref.JT 6 |) is obtained from the neutron /3-decay. 
We get 03(4 GeV 2 ) = cns(^ GeV 2 ) g A l§ = 1.11 and with this value we shall 
compare 03 calculated from our fits. Another combination of Vi(Q 2 ) and 
T”(Q 2 ) is equal to: 


T?(Q 2 ) + T™(Q 2 ) = 5a 8 (Q 2 )/18 + 2As(Q 2 )/3 


( 12 ) 
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with as = c NS (Q 2 )g8 and where g$ = 0.58 ± 0.03 |L6| is obtained from the 
hyperon /3-decays. Knowing cns{Q 2 ) we can calculate ds(4GeV 2 ) = 0.515 
and with this number we shall compare the results obtained from our fit. If 
we have had very precise experimental data and in the whole x range there 
would be no problems with determination of polarized quark distributions. 
Unfortunately that is not the case yet. Actually, from the experiment we have 
information on Tf and T”. The combination Uj'-T" is directly connected to 
qa experimental quantity modified by QCD corrections. On the other hand 
Ti+r^ is the combination of a% and As and it comes out that the fits are 
not sensitive enough to determine as and As separately in a stable way. The 
value of as and As are different for models I and II and different subsets of 
data. To stabilize the determination of parameters we assume in addition 
that as = 0.515 with 0.1 as artificial theoretical error. We will not present 
the results of all our fits. Some examples are given in the Table 1. 


Table 1 


The first moments of polarized distributions (see eqs.(9) and (10)). The 
strange sea polarization As is connected to the total sea polarization by the 
relation: As = 0.196AM. We present figures for two type of models: I and 
II. We have made our fits taking different spin asymmetries, namely: for 
proton (p) (P stands for all proton data with the inclusion of newest SLAG 
E143 points), neutron (n) and deuteron (d) target. 



i — i 

p 

3 

HH 

M 

w 

p 

II(p,d) 

H(p,n) 

H(P,n) 

r? 

0.156 

0.146 

0.145 

0.148 

0.144 

rn 

1 1 

-0.045 

-0.055 

-0.076 

-0.038 

-0.042 

9 a 

1.205 

1.201 

1.330 

1.114 

1.112 

9s 

0.508 

0.512 

0.515 

0.523 

0.478 

AS 

0.384 

0.297 

0.204 

0.373 

0.348 

AM 

-0.210 

-0.366 

-0.529 

-0.255 

-0.221 


It would be natural to use all the existing data for the determination of 
polarized quark distributions. We see that in this case fits I and II lead to 
similar results but formally fit I has one parameter more and higher y 2 per 
degree of freedom. From the Table 1 we see how integrated quantities: Tf, 
T", 03 , as, AS and AM depend on usage of different subsets of data. It seems 
that the closest value to the expected theoretical value for a 3 = 1.11 we get 


5 




















when we use the combination of all proton and neutron SLAC data. In the 
ratio as/as QCD corrections cancel so it can be expressed directly through the 
quantities known from low energy decay experiments i.e. g 8 / 9 a — 0.46±0.02. 
For the combination of proton and neutron SLAC data the ratio calculated 
from the fit ( a 8 /a 3 = 0.47) is closest to the value expected from experiment. 
We get for this fit y 2 = 20.47. If we add to this value the figure 6.93 
corresponding to deuteron contribution (from the SMC experiment) we get 
the value 27.4 what is very close to the the y 2 value obtained from the fit 
to all ( proton, neutron and deuteron) data, namely 27.1. If one starts, on 
the other hand, with the fit to proton and deuteron data and adds to the 
y 2 value the number corresponding to neutron contribution one gets y 2 = 
37.1. It seems that it is possible to get a satisfactory deuteron asymmetry 
using proton+neutron data and the result is not as good when we start 
with the combination of proton and deuteron data (in this case <23 = 1.33, 
the value which is far too big). So we have decided to use combination of 
proton and neutron ( 3 He) spin asymmetry data to calculate the parameters 
of the polarized quark distributions. We present in Figs.(la,lb,2 and 3) the 
comparison of our fit II(p,n) with the experimental asymmetries for proton 
(la,lb) neutron (2) and deuterium (3) target. In the last figure our curve 
should be treated as a prediction, because we do not take deuteron data into 
account in this fit. The parameters given in the first row of the Table 2 
correspond to the considered model, namely the one from the fifth column 
of Table 1. 


Table 2 


The coefficients of polarized distributions (see eqs.( 6 )) for type II models. 
The second row corresponds to the figures gotten in a fit with the inclusion 
of the new proton data from the E143 experiment at SLAC. 



CL\ 


a 3 

61 

h 

h 

c 2 

d 

H(p,n) 

0.924 

-3.237 

11.30 

-0.029 

-0.163 

-1.644 

-0.993 

-0.015 

H(P,n) 

0.929 

-3.005 

10.24 

-0.066 

-0.064 

-1.644 

-0.861 

-0.013 


The obtained quark distributions lead to the following integrated quanti¬ 
ties: A u = 0.91 (A u v = 1.04), Ad = —0.33 (A d v = —0.18) and As = —0.07. 
It gives the amount of sea polarization AM = —0.26. We see that our inte¬ 
grated quantities (T? and T”) differ slightly from the values quoted by the 
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experimental groups. The experimental figures are calculated directly from 
the experimental points with the assumption of Regge type behaviour at 
small x. On the other hand our polarized quark distributions satisfy all the 
constraints taken implicitly into account in fits to unpolarized data. 

We have also tried to consider gluonic contributions which were proposed 
as a solution of spin ” crisis” . We have taken into account gluons by consid¬ 
ering effective polarized quark distributions A q 6 ™ = A q — ^-AG where 

A G(x) = 1 — x ) 5 ' 3 , (13) 

with a new / constant which was fitted. In the basic fit (II(p,n)) we have 
got a s AG/2ir = —0.08 ( which corresponds to AG = —1.84 for a s = 0.28). 
This means that the sign of the gluonic contribution is opposite to what is 
expected theoretically and effective AT, 6 ” is bigger than AS coming from 
quark distributions. We conclude that from the point of view of fitting we 
do not need gluonic contribution. 

The Fig. (4) shows the comparison of g\(x) calculated from our fit with 
experimental points. We observe the substantial growth of g\ for small x 
values in our model. 

The last column of Table 1 corresponds to the inclusion of recent data 
for spin asymmetries from the experiment E143 at SLAG. We see that the 
inclusion of new data changes integrated quantities only slightly what we 
consider as a positive fact. In the case of the good fit the additional experi¬ 
mental data should not change parameters in a drastic way. The parameters 
corresponding to this fit are given in the second row of Table 2 whereas the 
fitted curve is presented in Fig.(5) together with the one for fit II(p,n). We 
also would like to point out that with the new data from E143 included the 
fit of the type I model coincides with that of type II, which means that the 
more singular behaviour of the sea is eliminated by fitting the formulas to 
the experimental data. However, one should be aware of the fact that in the 
last fit we have much more proton than neutron data points. 

Now, we would like to make some comments about the x —> 1 behaviour 
of valence quark distributions. Looking at the data points for proton spin 
asymmetry the value close to 1 at x ~ 1 is preferred, whereas for neutron and 
deuteron case values close to 0 seem to be natural (such observation is fragile 
due to the big experimental errors in this x region). In our approach we can 
give predictions for the behaviour of polarized quark distributions and spin 
asymmetries in the x —> 1 limit. In our fit we get Ai = A™ = Af = 0.89. 
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Starting from the new, improved version of the MRS fit [§ to the unpo¬ 
larized deep inelastic data we have made a fit to proton, neutron ( 3 He ) and 
deuteron spin asymmetries in order to obtain polarized quark parton distri¬ 
butions. We have discussed two kinds of models with more and less singular 
sea contribution at small x (both models being a possible consequence of the 
unpolarized distributions) and different combination of proton, neutron and 
deuteron spin asymmetries data. To stabilize the fits we add the experimen¬ 
tal information on octet quantity as- We have calculated the parameters of 
the polarized quark distributions from the combination of proton and neu¬ 
tron spin asymmetries data (without the SMC deuteron data). The inclusion 
of the new proton SLAC data from E143 experiment modifies our fit only 
slightly, which is a positive feature. We do not need gluonic contributions to 
be taken into account, i.e. the fit with gluons is worse. The new consistent 
data for spin asymmetries (mainly for deuteron and neutron) could help us 
to determine polarized quark distributions and to resolve the doubts about 
the x —> 0 and x —> 1 behaviour of such functions. 
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Figure captions 

Figure la The comparison of spin asymmetry on protons (data points are from 
SLAC (E80, E130) and CERN (EMC,SMC) experiments) with the 
curve gotten from our fit II(p,n) (eqs.(6,7) and the first row of Table 
2 ). 

Figure lb The same as in figure la but with x in logarithmic scale. 

Figure 2 The comparison of spin asymmetry on neutrons (SLAC E142 data) 
with the curve gotten from our fit. 

Figure 3 Our prediction for deuteron asymmetry compared with the SMC data. 

Figure 4 The data for g{{x) structure function with the curve gotten using the 
parameters of fit II(p,n). 

Figure 5 The same as in Fig.(lb) but with the new experimental points (E143) 
from SLAC included. The dashed curve corresponds to the fit lI(P,n), 
which takes into account the E143 experiment data. 
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